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Abstract 



Interconvertible enzyme cascades, exemplified by the mitogen activated protein kinase (MAPK) 
cascade, are a frequent mechanism in signal transduction pathways. There has been much 
speculation as to the role of these pathways, and how their structure is related to their 
function. A common conclusion is that the cascades serve to amplify biochemical signals 
so that a single bound ligand molecule might produce a multitude of second messengers. 
Some recent work has focused on a particular feature present in some MAPK pathways 
- a negative feedback loop which spans the length of the cascade. This is a feature that 
is shared by a man-made engineering device, the feedback amplifier. We propose a novel 
interpretation: that by wrapping a feedback loop around an amplifier, these cascades may 
be acting as biochemical feedback amplifiers which imparts i) increased robustness with 
respect to internal perturbations; ii) a linear graded response over an extended operating 
range; iii) insulation from external perturbation, resulting in functional modularization. We 
also report on the growing list of experimental evidence which supports a graded response 
of MAPK with respect to Epidermal Growth Factor. This evidence supports our hypothesis 
that in these circumstances MAPK cascade, may be acting as a feedback amplifier. 
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Introduction 



The ability of biological cells to receive and respond to signals, such as environmental conditions, 
mating potential or developmental cutes is considered a fundamental characteristic of life. In 
multicellular organisms, cells will often receive a flood of signals aimed at regulating their 
behavior. It is not surprising therefore to discover that cells have evolved highly elaborate 
and complex networks of 'signaling' pathways whose role is to coordinate and integrate this 
information and to elicit a suitable response. One of the most surprising results to have come 
to light in the study of these signaling pathways is the re markable degre e of evolutionary 
conservation that exists among widely different organisms (jMarshalll . 1 19941 ). The family of 
MAPK (mitogen- activated protein kinase) pathways, in particular, is highly conserved. 

The MAPK pathwa ys, besides being highl y conserved, are common components in signal 
transduction pathways ( Chang <fe Karinf . 2001 ). Virtually all eukaryotic cells that have been 
examined (ranging from yeast to man) possess multiple MAPK pathways each of which re- 
sponds to multiple inputs. In mammalian systems MAPK pathways are activated by a wide 
range of input signals including a va riety of growth factors and environmental stress es such 
as osmotic shock and ischemic injury ( Kvriakis Sz Avruch . 2002 ; Gomperts et al , 20021 ). Once 
the MAPK pathways have integrated these signals, they coordinately activate gene transcrip- 
tion with resulting changes in protein expression leading to cell cycling, cell death and cell 
differentiation. 

Modern methods of genetic manipulation h ave allowed ext ensive investigation of MAPK 
pathways in budding yeast. Through this work ([Marshall 1 19941 ) the notion of multiple parallel 
MAPK sig naling pathway s was developed. At least six S. cerevisiae MAPK cascades have been 
identified ( Banuettl . Il998l ). Many of these pathways receive multiple signals, some of which are 
shared with other MAPK cascades. This signalling architecture allows for cross-talk between 
the input stages of the MAPK pathways. The emerging picture is that of signal transduction 
pathways consisting of least two functional layers: a decision- making layer responsible for re- 
ceiving and integrating the inputs, which are then channeled to one or more MAPK stages in a 
second layer for final delivery to the target, in this case transcriptional modulation (see Figure 
1). 
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Figure 1: MAPK Pathways 
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In this paper we address the functional nature of the MAPK pathway stage. What properties 
do MAPK pathways possess which might suggest their role in this biochemical computation? 
There are a number of characteristics which are common to all MAPK families and have a 
bearing on how we should answer this question. Firstly, the structure of these pathways can 
give rise to extremely high gains on the input/output response, of the order of hundreds of 
thousands in some cases. Secondly, MAPK cascades always appear as the final stage of a 
signal transduction pathway. To our knowledge there has never been reported any further 
networks between the last stage of a MAPK cascade and transcriptional modulation. Finally, 
these cascades appear t o have little interaction with the rest of the system beyond their basic 
input/output behavior ( Widmann et al . 19991 ). In fact, this functional separation may be 



achieved in part by a physical separation, as some evidence suggests that the components of 
MAPK cascades are held together by scaffolding p roteins which may provide a spatial separation 
from the rest of the signal transduction network ( Widmann et al. . 19991 ; Garrington Sz Johnson . 
19991 ). 

These observations lead to a satisfying interpretation of MAPK cascades as modular com- 
ponents of signal transduction pathways which act to amplify thei r inputs. Others have also 
noted the possibility that MAPK cascades can act as switch de vices (IHuang &: Ferrelll . Il99fih . or 
even delay systems ([Nelson et all 120041 ; llhekwaba et all 120051 ) . In this paper we wish to inves- 
tigate the hypothesis that MAPK pathways, may in some circumstances be acting as feedback 
amplifiers. 

The paper is organized as follows. The structure and function of MAPK cascades is con- 
sidered first. The concept of Feedback amplifiers is introduced next, where the analysis is kept 
simple by considering linear systems. The benefits of feedback in the MAPK cascade are illus- 
trated in the next section, both by simulation results and by analysis of the linearization of the 
model. A discussion of experimental evidence is included in the last section. 
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Figure 2: Schematic of a Generic MAPK Pathway 
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Structure and Properties of Cascades As discussed in the introduction, MAPK cascades 
are common components of extracellular signal-regulated pathways. Experimental r esults sug- 



gest that these cascades are highly conserved across a wide spectrum of organisms ([Marshall 



1994 ). In addition, the basic structure of the pathway in each MAPK family is essentially the 



same. This implies that the functional roles played by each family is the same, or very similar. 
The architecture involves three protein kinases (see Figure 2) , each of which activates the kinase 
below only after being phosphorylated from above. 

A compell ing explanation of the function of these cascades is that they provide a switch- 
like response (jChock fc Stadtmanl . Il977l : iGoldbeter fc Koshlandl . Il984l ) which may be used in 
decision making, for example in development. It has been appreciated for many years that 
even a single cycle can elicit switch- like behavior under condition s in which the cycle reac- 



tions are operating near saturation (Goldb eter &: Koshlandl. 11981 ). A more general analy- 



sis was made by Small and Fell (IFell &; Smalll Il98fil : ISmall fc Fell . Il99(ll ) and Cardenas and 
Cornish-Bowden ( Cardenas Cornish-Bowden , 19891) using metabolic c ontrol analysis and sub- 
sequently extended by Ka hn and Westerhoff ( Kahn Sz Westerhofi . 1 99 ll ) and Kholodenko et al. 

19971 ). Arguments for switch- like behavior are supported by experimental 



( Kholodenko et al 



evid ence that some MAPK cascades may exhibit positive feedbacks leading to bistable behav- 
ior. (|Ferrelj liooi : iBhalla et ali l2002h 

During a switching action, the 'output' of the system (the concentration of activated kinase 
at the end of the cascade MAPK-PP) is switched from its basal (non-stimulated) le yel to satura- 
tion (i.e. 100% activation), as the 'input' (e.g. ligand level) crosses a threshold value ( Huang Ferreli 
1996l i This behavior is often described by the analogy to a cooperative enzyme - the cascade 
acts as a single enzyme with a large Hill coefficient. Along with its function as a memoryless 
switch, the cascade also amplifies the input signal, since the pool of MAPK typically has a 
concentration orders of magnitude higher than the ligand levels which activate the response. 

In addition studies have i ndicated the prese nce of one or more negative feedback loops sur- 
rounding the MAPK cascade ( Downward . 19961 ). This feedback is mediated by phosphorylation 
of SOS by ERKPP which in turn results in dissociation of the ShcGS complex thus interfering 
with Ras activation of the ERK pathway. 

In this paper we suggest a novel interpretation of the role of negative feedback in a MAPK 
cascade, inspired by analogous situations in electrical engineering. Amplifiers are some of the 
most ubiquitous components in electrical circuits. All amplifiers necessarily saturate above and 
below certain input levels. They are designed so that their 'active range' covers the levels of 
input of interest. When this active range is small, the result is switch-like behavior, suggesting 
an analogy with MAPK pathways. 

While amplification of signals is a common function of electrical circuits, amplifiers are 
rarely used in isolation. Rather, negative feedback is introduced from output to input, resulting 
in a feedback amplifier. These devices represent the key component of all (electrical) analog 
computation, since they allow the construction of functional operational amplifiers. The benefits 
of wrapping feedback around an amplifier may not be intuitively obvious, and indeed came as 
a surprise to many people when the idea was introduced. The primary improvements achieved 
with feedback are i) increased robustness with respect to internal perturbations; ii) a linear 
graded response over an extended operating range; iii) insulation from external perturbation, 
resulting in functional modularization. 

In the electronics industry these properties are exploited in the manufacture of amplifiers. 
Cheap, low tolerance, high-gain amplifier components are mass produced. These amplifiers are 
not suitable for use alone, but rather are coupled with high tolerance passive resistor components 
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which implement feedback and thus significantly improve the characteristics of the amplifier. 

We present the hypothesis that evolution may have hit upon the benefits of the feedback 
amplifier, and is implementing it in the MAPK cascade. In doing so we build on the work 
of others who have su ggested that negat i ve feedback may serve to improve the functi on of 
biochemical amplifiers ( Bhalla et al . 20021 ; Cinquin Sz Demongeot . 2002 ; Savageau . 119761 s ! . 

In its active range, the cascade would act solely as an amplifier - providing a faithful am- 
plification of the input presented to it from the mechanisms upstream. To these upstream 
mechanisms (which we do not address) would then be relegated the all-important task of in- 
tegrating the myriad extracellular signals received by the cell. The MAPK cascade then plays 
the role of amplifying the results of those biochemical calculations to the point where they 
can produce a response in the cell. This is an obvious strategy, again exemplified in electrical 
engineering. Analog computers perform their calculations on tiny currents and voltages then, 
when finished, amplify the results so that they may be useful. The same energy-conserving 
principle may be used in the cell. Signals may be processed at small concentrations (with cor- 
respondingly low needs for biosynthesis of enzymes). Only after the computation is complete 
and the appropriate response has been determined is it necessary to amplify concentrations to 
the point where they will affect the activity of the cell. 

That the MAPK pathway serves as the tail-end for wide variety of different signaling path- 
ways indicates that it may be an example of modular design. This design principle produces 
component which can play their role in a number of different situations. It may be that the 
MAPK cascade serves as a "plug-and-play" amplifier with specificity determined by inclusion 
of the appropriate "adapter", serving to connect the cascade to its input and output. This 
interpretation is supported by the finding that the "internal" mechanism of the MAPK cascade 
(i.e. the phosphorylation of MAPKK and MAPK) are highly conserved while the "connectors" 
(phosph orylation of MAPKKK and the target o f MAP K-PP) vary according to the role of the 
cascade ( Marshall . 1994 ; Garrington fc Johnson . 19991 ) 

This conservation of structure and mechanism is a great boon to the modeler. Based on 
this evidence we can have some confidence that analysis of a model of the cascade may provide 
insight into the biochemistry of a wide variety of signal transduction pathways. With that hope 
in mind, we present a n analysis of a model o f the Raf/Ras/MEK/E RK pathway presented 
m dKholodenkol . l2000h and based on work in (jHuang &; Ferrell Il99fj ). The model does not 
include the important mechanisms which occur upstream from the cascade. 

Despite the evidence for conservation of the cascade, there are variations in structure and 
mechanism which may or may not correspond to differences in pathway response. Common vari- 
ants are cascades which are longer or shorter than three kinases, or in which the kinases are phos- 
phorylated at only one site. Previous work s uggests that such d ifferences may not have much 
of an effect on the function of the pathway ( Brightman Fell . 2000l : Goldbeter Sz Koshland . 
19841 ). A more extreme variation is a cascade in which the activation is not performed by 



kinases, but rather by methylation or acetylation. There are fewer experimental results in such 
cases, but work by modelers suggests that an analogous functionality would be expected. In 
particular, in this paper we are concerned primarily with the cascade's qualitative role as an 
amplifier, which will likely be maintained over a wide variety of biochemical implementations. 



Feedback Amplifiers Rather than species concentrations, the amplifiers used by electrical 
engineers are designed to amplify current or voltage in an electrical circuit. A simple example 
is the voltage amplifier (i.e. voltage controlled voltage source) which samples the voltage in one 
part of a circuit and produces a proportional voltage in another. 
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Figure 3: Linear Feedback System: d = disturbance; 
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Amplification is one of the most fundamental tasks one can demand of an electrical circuit. 
One of the challenges facing engineers in the 1920's was how to design amplifiers whose per- 
formance was robust with respect to the internal parameters of the system and which could 
overcome inherent nonlinearities of their implementation. This problem was especially critical 
to the effort to implement long distance telephone lines across the U.S.A. 

These difficulties were ov ercome b y the introduction of the feedback amplifier, designed 
in 1927 by Harold S. Black ( Mindell . 2000l ). who was an engineer for Western Electric (the 
forerunner of Bell Labs). The basic idea was to introduce a negative feedback loop from the 
output of the amplifier to its input. At first sight, the addition of negative feedback to an 
amplifier might seem counterproductive. Indeed Black had to contend with just such opinions 
when introducing the concept - his director at Western Electric dissuaded him from following up 
on the idea, and his patent applications were at first dismissed. In his own words, "our patent 
appli cation was treated in the same manner as one for a perpetual motion machine" (jBlackl . 
1977h . 

While Black's detractors were correct in insisting that the negative feedback would reduce 
the gain of the amplifier, they failed to appreciate his key insight - that the reduction in 
gain is accompanied by increased robustness of the amplifier and improved fidelity of signal 
transfer. This trade-off between gain and system performance can be elegantly demonstrated 
by considering linear systems, to which we now turn. 



Results 

Consider the block diagram in Figure [3j We will consider only the steady-state behavior of the 
system, so we take the input u, the output y, the error e, and the disturbance d to be constants. 

*As Horowitz and Hill put it in (|Horowitz fc Winfieldl . Il~990l ) 'Negative feedback is the process of coupling the 
output back in such a way as to cancel some of the input. You might think that this would only have the effect 
of reducing the amplifier's gain and would be a pretty stupid thing to do.' 
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To begin, consider the case of no disturbance (d = 0). Assume that both the amplifier A and 
the feedback F act by multiplication. Without feedback (i.e. with F = 0), the system behavior 
is described by y = Au, which is an amplifier with (open loop) gain A. 

Introducing feedback, the behavior of the system is as follows. Prom the diagram 

y = Ae and e = u — Fy. 

Eliminating e, we find 

Au ID 



a l + AF 

Calling G = the system (or closed loop) gain, we have simply y = Gu. Comparing G with 

A, it is immediate that the feedback does indeed reduce the gain of the amplifier. Further, if 
the loop gain AF is large (AF S> 1), then 

r~ A 1 
~ ~AF ~ F' 

That is, as the gain AF increases, the system behavior becomes more dependent on the feedback 
loop and less dependent on the amplifier itself. We next indicate three specific consequences of 
this key insight. 

Sensitivity to internal parameter variation Since the closed loop gain G depends on the 
amplifier gain A, a small change AA in the gain A will produce a change AG in G. Considering 
the ratio of relative changes we find 

dG A 1 AG ^ 1 AA 

dAG ~ l + AF S ° ~G ~ l + AF A ^ 

provided AA is small. That is, a 1% change in the open loop gain A leads to a roughly 
pr^p % change in the system gain. We see that as the loop gain AF is increased the system 
becomes less sensitive to perturbations in its internal structure. In electrical circuit design this 
might alleviate the problems caused by temperature fluctuations or aging components. In a 
biochemical network the result will be a system whose function is not sensitive to fluctuations 
in enzyme concentrations. It should be noted that this reduction in sensitivity comes at a 
price - as the loop gain is increased the relative sensitivity to changes in F grows to unity. 
As mentioned, this performance constraint can be met by coupling low tolerance amplifiers to 
carefully tuned feedback elements. At this point we can only speculate as to whether nature 
has hit upon the same design strategy 

Sensitivity to disturbances in the output Suppose now that a constant disturbance d 
affects the output as in Figure [3l The system behavior is then described by 

y = Ae — d where e = u — Fy. 

Eliminating e, we find 

Au - d 
V ~ l + AF' 



S 



The sensitivity of the output to the disturbance is then 

dy = 1 
dd~ 1 + AF' 

Again, we see that the sensitivity decreases as the loop gain AF is increased. Such a disturbance 
could be caused by an increased load in a downstream part of an electrical circuit or removal 
of end product from a biochemical pathway. 

This property is of particular interest because the last stage of MAPK, that is MAPK-PP 
(see Figure 2) has to migrate to the nucleus in order to elicit a response. This diffusion is 
effectively a load on the MAPK circuit. Without feedback such a load would have a deleterious 
effect on the functioning of the MAPK pathway. As long as there is a pool of unphosphorylated 
MAPK, the feedback is able to compensate for this increased load. It is only when that pool 
dries up (as the amplifier is saturated) that the feedback is unable to provide any benefit. This 
property also has the natural benefit of automatically modularizing the network into, effectively, 
a single functional unit. 



Improved linearity of response over extended operating range Consider now the 
case where the response of the amplifier is a nonlinear function of the input signal, so the 
open loop response is y = A(u), e.g. y = 10^j). Taking the feedback again to act simply by 
multiplication, the behavior of the system G (now also a function of u) is described by 

y = G{u) = A(e) where e = u — Fy = u — FG(u). 

Differentiating we find 

de 

G\u) = A'(e)=- = A'(e)(l - FG'{u)). 
du 

Solving for G'(u) we find 

A'(e) 



G'{u) 



l + A'(e)F 

We find then, that if A'(e)F is large (A'(e)F > 1), then 

G'(u) » 1 

so, in particular, G is approximately linear, as its derivative is approximately constant. In 
this case, the linear feedback compensates for the nonlinearities in the amplifier A(-). Another 
feature of this analysis is that the slope of G(-) is less than that of A(-), i.e. the response is 
"stretched out". For instance, if A(-) is saturated by inputs above and below a certain "active 
range", then G(-) will exhibit the same saturation, but with a broader active range. 

A natural objection to the implementation of feedback as described above is that the system 
sensitivity is not actually reduced, but rather is shifted so that the response is more sensitive 
to the feedback F and less sensitive to the amplifier A. However, in each of the cases described 
above, we see that it is the nature of the loop gain AF (and not just the feedback F) which 
determines the extent to which the feedback affects the nature of the system. This suggests 
an obvious strategy. By designing a system which has a small "clean" feedback gain and a 
large "sloppy" amplifier, one ensures that the loop gain is large and the behavior of the system 
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is satisfactory. Engineers employ precisely this strategy in the design of electrical feedback 
amplifiers, regularly making use of amplifiers with gains several orders of magnitude larger 
than the feedback gain (and the gain of the resulting system) . 

In summary a feedback amplifier provides the following desirable characteristics: 

1. Increased robustness with respect to internal perturbations. 

2. Insulation from external perturbation, resulting in functional modularization. 

3. A linear graded response over an extended operating range. 

It is hard to overstate the importance of the feedback amplifier in circuit design. One of 
the most important components in analog electronics is the operational amplifier, or op-amp. 
The original concept of the op-amp came from the field of analog computers in the 1940s 
where extremely high-gain DC amplifiers were employed to carry out analog computations 
such as addition, multiplication, differentiation and most important of all, integration. The 
operating characteristics of a particular op-amp were determined by the feedback elements. By 
changing the arrangements and types of feedback elements different analog operations could be 
implemented. Thus the same amplifier was able to perform a variety of operations with only 
small changes to the feedback elements. 

One of the key characteristics of an op-amp is the extremely high gain between the input 
and output signals. By analogy we can see that the basic MAPK cascade pathway, without the 
feedback loop, has the same high gain characteristic. In the next section we shall illustrate how 
a high gain MAPK cascade with negative feedback will act in an exactly analogous manner to 
the feedback amplifier. 



Interconvertible Enzyme Cascades as Feedback Amplifiers Having described the util- 
ity of negative feedback in electrical amplifiers, we now demonstrate that the same benefits 
can be reaped when feedback is wrapped around a biochemical ampli fier such as the MAPK 
cascade. Ou r anal ysis is based on a model from Ferrell and Kholedenko (jHuang &; Ferrelifl99r^ : 
Kholodenkd . |2000| ) . wit h a negative feedback of the form described in previous work incorporat- 
ing si milar mechanisms ( Asthagiri Lauffenburger , 2001 ; Kholodenkd . 200C ; Brightman Sz Fel3 . 



2000). Note that our conclusions do not depend explicitly on this particular model but instead 



is a function of the high gain achieved by the cascade of cycles and the presence of the negative 
feedback. The model simply serves as a means to illustrate the concept. The network is shown 
in Figure [2j We begin by describing some simulations which show that this nonlinear system 
exhibits behavior similar to that described above for feedback amplifiers. 

At the basal level of no feedback, the steady state behavior of the pathway is shown 
in Figure 01 As discussed above, the cascade acts as a switch - the concentration of acti- 
vated MAPK climbs quickly from its basal le vel to saturation (i.e. complete activation) at a 
parti cular "threshold" concentration of signal (jGoldbeter <fe Koshlandl . Il984l ; iHuang &: Ferrell 
199fih . For inputs near the threshold level, the pathway acts as an amplifier, with the out- 



put level of activated MAPK far exceeding the signal concentration. (This is a feature which 
is strong ly dependent on the increasing concentrations of kinase s at s uccessive levels of the 
cascade ( Goldbeter Koshland . 19841 ; Asthagiri fc Lauffenburger , 2001 ).) 



Figure [5] shows the input/output (Sig to MAP K-PP) behavior of th e system as the feedback 
gain (1/Ki) is increased. Here we see, as noted in (jBhalla et all hooi ) the 'S'-shaped response 
curve is stretched and straightened. As the feedback is increased, the operating range of the 
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Figure 5: Response for various levels of feedback gain 
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Figure 6: Effect on response on a change in V\ (See equation (4)). In the first panel, pairs 
of graphs are shown, each pair represents a wild type and a perturbation to V\. Individual 
pairs are measured under differing feedback strengths (Solid: no feedback (Ki = 10000), Dash: 
moderate level of feedback (Ki = 600), Dash-Dot:strong feedback (Ki = 200)). The second 
panel shows the difference in response between each pair. Note that as the feedback increases, 
the difference becomes more linear. 



amplifier is increased and the response is more linear over that range. It is also clear that 
the overall gain is decreasing, since the resp onse for each parti cular input level diminishes as 
the feedback gain increases. It is argued in dBhalla et all l2002h that this behavior allows the 
cascades to be bi-functional, acting as a switch when the feedback is low and as an amplifier 
when the feedback is high. From the point of view of the amplifier characteristics, a switch can 
be thought of as an amplifier with very short operating range. An increase in the feedback gain 
increases the range over which the cascade acts as an amplifier. 

Figure [6] shows the effect of parameter variation on the system at various levels of feedback 
intensity. The response of the system is plotted against signal level for three different levels 
of feedback gain. In each case, the response is also shown after a perturbation in the kinetic 
parameter V\ (a 20% decrease). Also shown is a plot of the difference in response between the 
nominal and perturbed systems. The feedback is again having a "stretching" effect. While the 
"total" error is not appreciably diminished, the effect of the feedback is to spread that error 
out across the increased operating range of the amplifier. The result is a significant decrease 
in the change in response for each particular input level within the active range. Similar result 
hold for perturbations in the other model parameters. 



Linearized model In addition to simulations of the complete model, analysis of a lineariza- 
tion is also useful in drawing analogy to the linear feedback amplifier described above. The 
model was linearized about a nominal operating point (steady state) with moderate feedback 
gain (Kj = 600) and an input level in the operating range of the amplifier (Sig = 3). The 
linearization takes the form 

x = Ax + Bu - F(Ki)x 
where the state variables x and the input u describe the deviation of the enzyme concen- 
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trations and signal level from their nominal values. This is a slightly different form than the 
linear amplifier discussed above since the feedback acts on the system rather than directly on 
the input (i.e. the feedback influences the effect of Sig on MAPKKK rather than having a 
direct effect on Sig itself). In this case we expect the closed loop gain to have the form 

Aiu 

y 



1 + AoF 



where A\ is the gain from input to MAPKKK-P, and A2 is the gain from MAPKKK-P to 
MAPK-PP. In this case the feedback is described as 

4.87 



so that 



n*i) - (1 + i64 )2Ki 



Am A K. 2 + 328Kj + 26900 

Aiu —* (3) 



l + A 2 F(Ki) K 2 + (328 + 4.87^ 2 )Ki + 26900 

The linear (steady state) relationship between input and output (i.e. signal level and MAPK- 
PP) is given by 

[MAPK-PP] - [MAPK-PP] = ([Sig] - [Sig] )2180 + 329K, + 27000 
1 J L J0/ K 2 + 15600Kj + 27000 

where [MAPK-PP] and [Sig] are the nominal levels. This matches ([3]) with A\ = 2180, 
A 2 = 3140. 

Another way of characterizing the linear system is through its frequency response. Graphs 
of the response versus frequency (i.e. magnitude Bode plots) are shown in Figure [7] for various 
levels of feedback. This is the familiar response of a feedback amplifier - as the feedback 
increases, the gain of the system decreases, but the range of frequencies over which that gain 
is maintained is increased. This highlights another feature of adding negative feedback to 
an amplifier: an increase not just in the range of constant signals that are amplified, but a 
corresponding increase in the range of sinusoidal signals that are amplified. Although this 
amplification of oscillatory signals has not been seen as biologically significant, one cannot rule 
out that some cells may be using the MAPK cascade for this purpose, especially in light of the 
number of oscillatory signals which have been identified. 



Computational Model The model network is shown in Figure [2j The eight species concen- 
trations of interest are MAPKKK, MAPKKK-P, MAPKK, MAPKK-P, MAPKK-PP, MAPK, 
MAPK-P, MAPK-PP. Moiety conservation allows us to treat three of these as dependent vari- 
ables through 

[MAPKKK] = Ti - [MAPKKK-P] 
[MAPKK] = T 2 - [MAPKK-P] - [MAPKK-PP] 
[MAPK] = T 3 - [MAPK-P] - [MAPK-PP]. 

The reactions are described by Michaelis-Menten rate laws: 
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Figure 7: Frequency response of linearized model for various feedback gains 



(Vi + [Sig])[MAPKKK] 

(Kmi + [MAPKKK]) (\ + [ MA ^- pp ] \ 

V 2 [MAPKKK-P] 
Km 2 + [MAPKKK-P] 

V 3 [MAPKK] [MAPKKK-P] 
Km 3 + [MAPKK] 
V 4 [MAPKK-P] 
Km 6 + [MAPKK-P] 
V 5 [MAPKK-P] [MAPKKK-P] 
Km 5 + [MAPKK-P] 
V 6 [MAPKK-PP] 
Km 4 + [MAPKK-PP] 

V 7 [MAPK] [MAPKK-PP] 
Km 7 + [MAPK] 
V 8 [MAPK-P] 
Kmio + [MAPK-P] 
V 9 [MAPK-P] [MAPKK-PP] 
Km 8 + [MAPK-P] 
Vio [MAPK-PP] 
Km 9 + [MAPK-PP] 
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Nominal parameter values are chosen as 



Kmi = 10 Km 2 = 8 Km 3 = 15 Km 4 = 15 Km 5 = 15 
Kmg = 15 K1117 = 15 Kins = 15 Knig = 15 Kmio = 15 
Vi = 2.5 V 2 = 5.1 V 3 = 0.025 V 4 = 0.75 V 5 = 0.025 
V 6 = 0.75 V 7 = 0.025 V 8 = 0.5 V 9 = 0.025 Vi = 0.5 
Ti = 100 T 2 = 300 T 3 = 300 

The units for concentration and Michaelis constants are given in nM, maximum rates are 
expressed as nM and rate constants in s _1 . To investigate the behavior of the system under 
a load on the output, sequestration of MAPKK-PPin the nucleus is included in the model by 
the inclusion of an additional state variable nucMAPK-PP, described by 

[MAPK] = T 3 - [MAPK-P] - [MAPK-PP] - nucMAPK-PP 



Vll = fen [MAPK-PP] V12 = fei 2 [nucMAPK-PP] 
fen = 0.003 fcu = 0.01 

The model was linearized about a nominal steady state with moderate feedback level (Kj = 
600) and a signal level in the amplifier range (Sig = 3). The nominal values of the state variables 
are 

[MAPKKK-P] = 23.8 [MAPKK-P] = 43.9 [MAPKK-PP] = 21.7 
[MAPK-P] = 82.3 [MAPK-PP] = 164 

The linearized model takes the form 

x = Ax + Bu - F(Ki)x 



where 



[MAPKKK-P] \ 
[MAPKK-P] 
[MAPKK-PP] 



A 



/ [MAPKKK-P 
[MAPKK-P 
x = [MAPKK-PP 

[MAPK-P] - [MAPK-P] 
\ [MAPK-PP] - [MAPK-PP] / 

/ -0.0462 
0.00487 -0.00596 0.00822 
0.0186 0.00257 -0.00837 





u = [Sig] - [Sig] 








-0.00163 
0.0211 







-0.00339 
0.000859 



\ 




-0.00151 
-0.000233 / 



B 



( 0.694 \ 
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Discussion 

In this paper we have discussed an alternative hypothesis on the functional role of the MAPK 
pathway. In particular, we have discussed the possibility that MAPK is acting as a classic 
feedback amplifier so beloved of engineers. 



Experimental Evidence Although the previous theoretical and computational arguments 
offer a novel view of how the MAPK pathway may be operating, the popularly held view is 
that MAPK pathways are in fact operating as switch devices rather than generating a graded 
response as we have suggested here. Indeed in some cases a switch response makes a very 
satisfactory explanation, for example, oocyte m aturation presumably requires an all-or-none 
response to progesterone ( Huang Sz Ferrelll . 1996). However in other cases, t he situation is not 
so easily decided. Hazzaline and Mahadevan ( Hazzalin Mahadevan . 20021 ) have discussed at 
length the case for a graded response over a switch like response. They provide evidence for 
the graded regulation of immediate-early genes in cell culture which are controlled via MAPK 
pathways. Furthermore, the re is now clear evidence to supp ort graded responses in the yeast 



20011 ; lFerrelil2002l ) which is again based on a MAPK 



mating pheromone pathway (jPoritz et al 
cascade. 

There are also a growing number of papers which report a graded response of ERK (Mitogen 
activat ed protein kinase) w ith respect to EGF (Epidermal Growth Factor) stimulation. The first 
paper (jBhalla et aUl2002l ) describes a cell population study where the response was investigated 
using Platelet-Derived Growth Factor (PDGF), they observed a proportional response over 
a 10-fold concentration range of PDGF. This effect was explained in terms of the action of 
the negative feedback modulator, MKP (MAPK Phosphatase). Although this feedback only 
operates on one level of the three cascade structure, it still induces the property of linearization 
due to negative feedback as described in this paper. 

Much more convincing however is recent work carried out on single cell studies. It is 
becoming increasingly apparent that population level studies give only a rough indication of 
cellular dynamics, this is due to heterogeneity in the population. A stricking ex ample of this 



effect are the oscillations recently observed in p53 in response to DNA damage (lLahav et al. 



20041 ) . In population studies the oscillations appear damped (due to each oscillator being out 
of phase) , while in single cell studies the oscillations are of fixed ampli t ude. It is therefore of 
interest to note two papers (IWhitehurst et all . |2004j; iMackeigan et all l2005h which focus on 
single cell studies. In both cases they o bserved a clear graded r esponse of ERK stimulation 
with respect to EGF. In the first paper dWhitehurst et all liooi ) the concentration of dually 
phosphorylated ERK in individual HeLa cells was correlated with activation as a result of EGF 
and PMA stimulation. In both cases a clear linear response was observed. The second paper 
(jMackeiean et all , hood ) investigated the response of Swiss 3T3 fibroblasts to EGF and PDGF. 
The work used a FACS (Fluorescence Activated Cell Sorting) based approach to investigate the 
distribution of activation in the cell population. If the cascade operated in an ultrasensitivity 
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mode then the authors expected to see two cell populations corresponding to an off an on state. 
Instead the experiments showed intermediate activation of ERK in response to EGF. 

Although these experiments indicate that ERK responds linearly to stimulations they do 
not suggest how the linearity originates. We postulate in this paper that such a response 
would originate from ultrasensitivity coupled with negative feedback. To determine whether 
this is the case further studies would need to be undertaken, one possibility would be to remove 
the feedback loop and to remeasure the response. Further evidence of the role of the negative 
feedback could be obtained by investigating the robustness of the response to changes in proteins 
levels inside the loop. If the feedback loop is operating then such changes should have little 
effect on the network response. While such experiments would not prove conclusively that the 
cascade is acting as a feedback amplifier, they would certainly lend weight to that alternative. 



Netw ork Versatility The MAPK unit appears to be a very versatile unit (jBhalla et al. 



2002). Without any feedback, MA PK can act as a very high gain amplifier, s o high that 



is behaves in a switch like manner dHuang & Ferrell Il996l : iBagowski etal\ . \200i ). With the 



addition of feedback a number of new behaviors can emerge. With positive feedback, MAPK 
can act as a bistable switch, with negative feedback MAPK can act as a classic feedback 
amplifier. Moreover, with sufficient negative feedback the MAPK pathway can be made to 
oscillate dKholodenkol. l200Ch. Many of these behaviors have now been observed experiment ally 
(|Huang Ferrelj Il996l : bhalla et all 120021 : Iwhitehurst et all \2004 : iMackeigan et all I2005T ) . 

One of the great challenges in modern molecular biology (or systems biology) is to under- 
stand the functioning of large complex networks. Such networks are composed of thousands of 
reactions and without some way to simplify their description we will find them almost impossi- 
ble to understand. In electrical engineering, it is now almost routine to design complex circuits 
composed of millions of components, for example the latest Pentium 4 from Intel is reportedly 
composed of forty three million transistors. How, one may ask, do engineers cope with such 
enormous complexity? A large part of the answer like in modularization, that is, the network is 
broken down into smaller function units which have distinct and well defined behaviors. These 
in turn are modularized further, as appropriate, down to the level of the discrete transistor. 
Thus a hierarchy of function is described. By modularizing, each module becomes manage- 
able in terms of design and understanding. In reac tion networks, we ne ed to take the same 
approach. Although the concept of modularization (|Lauffenburgerl . 2000) in cellular networks 



— ~ 

(Milo et al. 


2002; 


— „ 

Ravasz et al.. 


2002; 


Lee et al.. 



20021 ). 



However we have one major problem which the engineers do not have. Whereas man-made 
devices are designed, natural systems are clearly not. Thus we need to be able to reverse engineer 
the modules. This is obviously no easy task, for example we do not even know whether natural 
selection would a utomatically generate a modul e base d structure. Experiments on artificially 
evolved systems ( Koza et al . 19991 : A fc Sauro . l2004h suggest that both modular and non- 
modular 'designs' can evolve, with modularity being the most common outcome. If this is the 
case for natural systems then we will be able to make progress. 

In this paper we describe the properties of one possible module, that is the MAPK module. 
Many other modules most probably exist. We are currently investigating other potentially 
functional units from both prokaryotic and euka ryotic systems. Others have already started 
using this appro ach (Smolen et al . 19981 . 2000ol Jfcl) and tw o excellent r e views in this area by 
Wolf and Arkin (|Wolf fc Arkinl . I2003T ) and Tyson et. al. (|Tvson et all I2003T ) are well worth 
reading. This area of research would probably gain much from looking at the sorts of devices 
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that electronic engineers employ, including devices such as switches, amplifiers, oscillators, 
frequency filters amplitude filters, noise filters and amplifie rs, combinator i al log ic, homeostats, 
rheostats, logic gates and memory elements (list taken from ( Wolf &: Arkinl . l2003l )). For Systems 
Biology to succeed, knowledge and experience from different disciplines clearly needs to come 
together and electrical engineering has probably a significant contribution to make in allowing 
us to understand complex biological networks. 



Materials and Methods 



All simulations were carried out using a combination of M atlab (http : //www . mathworks . com/) 
and the Systems Biology Workbench ( Sauro et al , 20031 ) . Graphs were generated using Matlab. 
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